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Abstract 

There are several new features in production, oscillations and detection of the 
atmospheric neutrinos of low energies, E ^ 100 MeV. The flavor ratio, r, of muon 
to electron neutrino fluxes is substantially smaller than 2 and decreases with energy, 
significant part of events is due to the decay of invisible muons at rest, etc. Oscilla- 
tions in two-layer medium (atmosphere - earth) should be taken into account. We 
derive analytical and semi-analytical expressions for the oscillation probabilities of 
these "sub-subGeV" neutrinos. The energy spectra of the e— like events in water 
cherenkov detectors are computed and dependence of the spectra on the 2-3 mixing 
angle, $23, the 1-3 mixing and CP-violation phase are studied. We find that varia- 
tions of 623 in the presently allowed region change the number of e— like events by 
about 15 — 20% as well as to distortion of the energy spectrum. The 1-3 mixing and 
CP-violation can lead to ~ 10% effects. Detailed study of the sub-subGeV neutrinos 
will be possible in future Megaton-scale detectors. 
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1 Introduction 



Studies of the atmosplieric neutrinos have been performed mainly at energies 
E > (0.1 — 0.2) GeV. There are only few results on neutrinos of lower energies: 
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E = (10— 100) MeV. The fluxes of these neutrinos [1] have been computed re- 
cently by several different groups [2,3], [4,5], [6], and there are some differences 
of the results of computations. 

The e— like events induced by the low energy atmospheric neutrinos have been 
detected by SuperKamiokande [7,8,9]. About 88 ± 12 events were produced 
by interactions of the atmospheric and directly, and 174 ± 16 events 
originated from decays of invisible muons. In turn, these muons are generated 
by the V/j^-Hux with typical energies (150 - 250) MeV [8]. LSD put only an 
upper bound on the z/g— flux: Fg < 5 ■ lO^cm^^s^^ [10] which is about 5 times 
higher than the expected value for the energy range 12 < £" < 26 MeV. 

As far as the oscillation effects are concerned, only the vacuum z/^ — z/^- oscil- 
lations have been taken into account in [7,8,9]. Also oscillations in two layer 
medium with constant densities relevant for low energy neutrinos have been 
considered [11]. 

The low energy atmospheric neutrinos were discussed, as a background for 
detection of the relic supernova neutrinos [12] as well as future SN neutrino 
bursts [13]. 

In this paper we perform detailed study of the oscillations of the "sub-sub- 
GeV" atmospheric neutrinos in the complete 3i/-context. The range below 
100 MeV offers rather rich oscillation phenomenology. Although presently the 
number of detected events is small, in future, new large scale experiments can 
accumulated large enough statistics to extract new interesting information. 

The paper is organized as follows. In sect. 2 we summarize properties of the 
neutrino ffuxes at low energies, as well as parameters of the neutrino trajec- 
tories. We present relevant oscillation probabilities inside the Earth in sect. 3. 
The probabilities of 3i/-oscillations in two layer medium (the atmosphere and 
the Earth) are derived in sect. 4. In sect. 5 we consider averaging and inte- 
gration of the probabilities over the angular variables. In sect. 6. we present 
the z/g, z/^ atmospheric neutrino fluxes at a detector. We then compute the 
numbers of e— like events as functions of energy, in water Cherenkov detec- 
tors induced by the direct i/g, Ue interactions (sect. 7) and via the invisible 
muon decays (sect. 8). We study dependence of observables on the oscilla- 
tion parameters. In sect. 9 we present the energy spectra of e— like events in 
the megaton-scale detectors. Discussion and conclusions follow in sect. 10. In 
Appendices A, B, C we present analytic and semi-analytic formulas for the 
oscillation probabilities. 
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2 Atmospheric v fluxes and trajectories 



Let us summarize properties of the neutrino fluxes with < 0.1 GeV. In flg. 1 
we show the produced (without oscillations) fluxes of the electron neutrinos, 
Fg , electron antineutrinos, Fg , muon neutrinos, F°, and muon antineutrinos, 
F^, as functions of the neutrino energy. The lines have been drawn using 
results of computations in rcf. [6]. The fluxes are averaged over the zenith and 
azimuthal angles. They depend on period of solar activity which influences the 
fluxes of primary cosmic rays. Shown in the flgure are the average values of 
neutrino fluxes during periods of the maximum and minimum of solar activity. 

The low energy neutrino fluxes are formed mainly in the decays of pions and 
muons in flight and at rest. The K-meson decays contribute less than 0.05%. 
The key energy scales (indicated by the dashed lines in the flgures) are the 
energy of muon neutrinos from the pion decay at rest, k, 30 MeV, and the 
end point of neutrino spectrum from the muon decay, E'^ = 53 MeV. 

Properties of the spectra can be summarized as follows. 

(1) For E > E^ = 53 MeV the i/g— and i^e— spectra are formed in muon 
decays in flight. 

(2) The bumps in i/g— and Pg— spectra below F^ = 53 MeV are due to the 
muon decay at rest. This contribution composes about 1/3 of the total 

flux at these energies. 

(3) The neutrino flux is slightly (5 — 10%) larger than the antineutrino flux: 

> Eg. The reason is that u^s originate from the chain of reactions 
TT"'" — > /I'^i'ij., /i'^ — > e'^i'gPi^, whereas Pg - from the conjugate reactions. 
Since the original cosmic rays are protons and nuclei, they overproduce tt"*" 
in comparison with 7r~, and consequently, the vr^ chain is more abundant. 

(4) For F > 53 MeV the muon (anti) neutrino spectra arc formed by the pion 
and muon decays in flight. Since both the n^— and tt^— decay chains 
produce equal number of and P^, the corresponding fluxes are approx- 
imately equal. The bump in the spectrum below E = 53 MeV with sharp 
edge at F = 53 MeV, is due to muon decay at rest. The peak at FJ = 30 
MeV originates from the pion decay at rest. Below 30 MeV the main 
contribution to the z/^— flux is from the muon decay, and about 38% of 
the flux is generated by the pion decay in flight with neutrinos emitted 
in non-forward directions. 

(5) Below Ej;^, the z/^^— flux is slightly larger than the u^— flux: F° = 1.05F°. 
The difference originates from the muon decay at rest and has the same 
reason as larger ffux of i>e'- comes from the chain tt /x"*" P^^. 

(6) According to [6] the tt— decay peak for P^ is larger than the f^— peak. 
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Fig. 1. The fluxes of z>^, z^^, and Ue neutrinos as functions of neutrino energy from 
[6]. The dashed vertical lines show the end point of the neutrino energy spectrum 
from the muon decay at rest, E^, and the neutrino energy from the tt— decay at 
rest, El. 

The flavor ratios 

'<^'®'> = iy(i:ej' ''^'''''■■> = Ffmrr 

play the key role in oscillations. Here is the zenith angle of neutrino tra- 
jectory. In what follow we will discuss the ratios averaged over the zenith and 
azimuthal angles. As follows from fig. 1, in the range E > E^^ the ratios equal 
r ^ 2.0 and f ~ 2.2. Below i?^ the ratios decrease with energy. For instance, 
in maximum, E ~ (35 — 40) MeV, one has r ^ 1.72 and f ^ 1.87; for E = 21 
MeV we obtain r ~ 1.65 and f ~ 1.77. This is important difference from the 
case of higher energy atmospheric neutrinos, where r ~ 2 and f ~ 2.2. 
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In fig. 2, we show the flavor ratio for the sum of neutrino and antineutrino 
fluxes: 

obtained from different computations. The ratio changes from 2 to almost 1.5, 
when the neutrino energy decreases from 100 to 10 MeV. Sharp jump of the 
ratio at 53 MeV is due to the muon decay spectrum which has maximum at 
the end point. 
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Fig. 2. The flavor ratio (2) obtained from results of computations of different groups: 
dashed line [3] , red line [6] , blue line [5] . The vertical lines indicate the neutrino en- 
ergy from the pion decay, (dashed), and the end-point of muon decay spectrum, 
EH (solid). 

There are differences between neutrino spectra presented by different groups. 
In fig. 3 we compare the sum of z/g and z/g- fluxes from three available com- 
putations. Notice that above 20 - 30 MeV the difference of spectra is about 
(10 — 15%) and the shapes of spectra are rather similar. To have an idea about 
the solar activity effect we show in fig. 3 the fluxes in maximum and minimum 
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of the activity computed by Bartol group [3]. The difference is (15 — 20)%. In 
what follows we will use the averaged over the solar cycle fluxes. 




0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 01 

E (GeV) 

Fig. 3. The (fe+Pe)- neutrino flux as function of energy for maximum (solid red) and 
for minimum (solid dashed) of the solar activity from [3] . Show are also the averaged 
over the solar activity cycle fluxes from [5] - Honda (green) and [6] - FLUKA (blue). 



In contrast to high energies, for low energy neutrinos one needs to take into 
account the oscillations driven by the 1-2 mixing in atmosphere. Indeed, the 
oscillation length is given by 

Amii V30MeVy ^ ^ 

which is comparable to the length of trajectory in the horizontal direction. 
The total length of neutrino trajectory from a production point to a detector. 
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L, is given by 



2R 



R cos^ Q^ + h- cos e^'^/i?^ cos2 6^ + 2i?/i + + 



2i? 



(4) 



where R is the radius of the Earth, /i ~ 20 km is the height in the atmosphere 
where neutrinos are produced, is the zenith angle. We neglect here a depth 
of detector below the surface of the Earth. 



Above the horizon, cos0,^ > 0, trajectories are in the atmosphere only. In the 
horizontal direction, cos Gj, = 0: L = \/2Rh + h"^ ~ y2Rh. For trajectories 
below the horizon cos O,^ < the length of trajectory inside the matter of the 
Earth equals 

= — 2i?cos©,^, (5) 
and for these trajectories the length in the atmosphere: 

L^(cos 0,^) = L(cos 0,^) — Lj„(cos 0,^) = L(— cos 0,^). (6) 

For trajectories not very close to the horizon we have ~ /i/ cos©:^ which 
corresponds to the flat atmosphere and L ^ 2^r{Rcos'^ + h). The length 
of trajectory in mater and in atmosphere become comparable. La = ~ 500 
km at I cos 0^1 ^Jh/2R = 0.04. 

The half-phases of oscillations in the atmosphere and the Earth equal 

= 7r— , =7r— , (7) 

''V I'm 

with ly and being the oscillation lengths in vacuum and in matter corre- 
spondingly. Notice that at low energies the phase in the atmosphere can not be 
neglected even for trajectories not very close to the horizon: for instance, for 
£; = 30 MeV and cos0^ = -0.3 we obtain sin0 = 0.2. For = 60 MeV, we 
have sine/) = 0.1. So, in general at low energies the vacuum oscillation phase 
can not be treated as a small parameter. The oscillations in atmosphere can 
be neglected for the muon neutrinos producing the invisible muons. 



3 Oscillation effects inside the Ecirth 



Let us first consider the main features of oscillations at low energies inside the 
Earth. The resonance energy equals 

Ami, cos 2^12 Qfi A M.v ( ^^'i ^ /2X)g/cm3\ /cos2^\ 
2K ^^^■^^'^1,7.3-10-W^J l^^J lvT42rj- 

(8) 
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So, the energy interval E < 0.1 GeV is at and below the 1-2 resonance. It is 
much below the 1-3 resonance: E/E^^^'^ < 0.015. As a consequence, 

• Matter effect on the 1-3 mixing is very small and can be neglected in the 
first approximation. 

• The third eigenstate approximately coincides with the third mass eigenstate 

• The state 1^3 decouples from dynamics and evolves independently. The two 
other states form 2i/— mixing system, and so the problem is reduced to 21/— 
problem. 

• Oscillations driven by the large mass split, Am|]^, arc averaged out in prob- 
ability due to integration over the angular variables and energy. 

In what follows we will quantify this picture and derive the relevant oscillation 
probabilities. Evolution of the neutrino flavor states, Vf = (ve, t^ix, '^t)^ ■, is 
described by the equation 

dvf (UM'^U^ A 

where the PMNS mixing matrix defined through Uf = U Vmass can be param- 
eterized as = U23VsUi3Ui2- Here = diag(l, l,e*''), V = diag{V, 0,0), and 
V — v^C^ne, and Uij is the rotation in the ij— plane onto the angle 

Consider new basis of states, u' = {j^e, ^'/x, ^3)^ , defined by 

Uf = U'u', (10) 

where 

U' = U23T5U13. (11) 
In this basis the Hamiltonian becomes 

H' = ^U,2M'U12 + UlVUl. (12) 

Taking into account that VS13C13 <^ Amli/2E we can perform the block- 
diagonalization of H' in (12) which leads to 

H' « ^UuM'uU + diag{Vci„ 0, Vsj,). (13) 

We use notations C13 = cos ^13, C12 = cos ^12, S12 = sin ^12, etc.. Very small 
term —V cl^sl^{2V E / Aml^) in the 11-element is neglected. (In this case the 
diagonalization is reduced to just omitting the off-diagonal elements in the 
second term of eq. (12).) Essentially, this approximation corresponds to ne- 
glecting the matter effect on 1-3 mixing. The matter correction to the 1-3 
mixing equals 9i3{2VE/ Aml^), and indeed, can be neglected [15]. According 
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to (13), the state z/3 dccoTiplcs from the rest of the system and evolves inde- 
pendently. The two other states form usual 2z/— system with the Hamiltonian 
H2 = i?2(Am|i, 6*12, Vc\'^). Therefore the S-matrix (the matrix of amplitudes) 
in the basis v' has the following form : 



S' 



A 



33 



(14) 



where 



A33 = exp(-i2(/)g^) , 20^ = 



2£; 



(15) 



and is the total distance traveled by neutrinos in matter. The amplitude 
describes transition v'^ — > v'^. The off-diagonal elements of S' equal each 
other: 

= K,, (16) 

as a consequence of the T-symmetry (the CP-conservation in this basis and 
symmetry of the Earth matter profile). The amplitudes A'^^, A'^^ and A'^^ arc 
solutions of the evolution equation with the Hamiltonian if2- As a consequence 
of unitarity, 

A' — A'* A'* — _ A' (-\7\ 

or ReAg^ = 0, that is, the transition amplitude is pure imaginary. 
According to (10) the -S"— matrix in the original flavor basis equals 

Sf = U'S'U'l (18) 
Therefore the — > i/p oscillation probability is given by 



P{u^^Uf,) = \{U'S'U'^)paf. 



(19) 



Using U' defined in (11) and S' from (14), and averaging the oscillations related 
to the third mass eigenstate we obtain explicitly 



2 

P(i/^ I/e) = C?3 -Si3S23e''^A'^e + C23A' + sj^clssl^, 



(20) 
(21) 



and for the inverse channel: Pl^e ^fj.) — P{^ix ^e){S S)- Finally, 

2 



Let us introduce three functions 



cl^A'„„ - Si3 cos S sin 26'23^eM + ^13^23^66 



^ '-13*23- 



(22) 



D = {P2,i?2,/2} 



(23) 
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(essentially the elements of the density matrix) as 

P2^\^^f = l-\AJ\ R,^Re{A:^AJ, h^\m{A:^AJ. (24) 

These functions satisfy the relation 

Pi + Rl + ll = P (25) 

which follows from unitarity of the S— matrix. The other properties of D— functions 
have been studied in [15]. 



According to (17) and (24) 



1 



Re(^';A,) = -^2, Re{A;^Ke) = -^{Ii-Ri)- (26) 

-'2 

Using these equalities and notations (24) the probabilities (20), (21) and (22) 
can be rewritten in terms of D as 



P(z/e-i^e)=4(l-P2) + 4, (27) 

P{,i^^i ^e) = cl^cl^P2 - sncls sin26'23(cos5i?2 + sin 5/2) 

+ sl,cl,sU2-P2), (28) 

(which coincides up to the sign of 5 with expression in the paper [15]) 0, and 



P(z/^ ^ z/^) = 1 - 0.5 sin^ 26*23 - C23P2 + 2si3 cos S sin 26'23(c23 - 513523)^2 + 



+ sfg sin^ 2^23 



2'5i3'S23 + 



+ .?34(2-P2). (29) 

The first line in this equation is a sum of the averaged standard 2z/-oscillation 
probability z/^ —>■ driven by the 2-3 mixing, the contribution from the 
1-2 oscillations for S13 = and the first (linear in S13) correction due to 
the 1-3 mixing (the "induced" interference [15]). The probability (29) is an 
even function of 5. Notice that in the second order the correlation of S13 and 
5: si3COS(5 is broken and these two parameters enter differently, which, in 
principle, opens a possibility to determine them independently. 

For antineutrinos we have the same expressions for the probabilities with sub- 
stitutions: 5 -6, P2 P2, R2 R2, h h, where P2 = P2{V -V), 
etc.. 



^ In our previous paper [15] we used different definition of amplitudes, A^^, for 
transition — > fe, instead of A^fj,. 
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Some insight into properties of the probabihties and their dependence on the 
osciUation parameters can be obtained using expressions in the constant den- 
sity case: 

= sin^ 2^^ sin^ (fT, = -- sin 4^™ sin^ (f)"", = - sin 2^^ sin 20"*. 

^ ^ (30) 

where subscript "c" refers to the case of constant density, and 0"* is the half- 
phase of neutrino oscillations in matter defined in (7). The probabilities for 
antineutrinos, P2, R2 and Jg, can be obtained from those in eq. (30) substi- 
tuting 9m Om and 0"* — > 0"*. The formulas for constant density (30) give 
good qualitative and in many cases quantitative description of the results. 

We compute the 2z/— probabilities R2 and I2 and similar probabilities for 
antineutrinos numerically. A very precise semi-analytical description of the 
probabilities can be obtained with the adiabatic Magnus expansion. In the 
first order of this expansion using the results of [16] (see eq. (78) and (76) in 
[16]) we obtain for the amplitudes in the v' basis 



A'gg = cos le cos 0"'' + i (cos le sin 0"'^ cos 2^°2 - sin h sin 291^ , (31) 

A'^^ = -i (sin le cos 29% + cos Ig sin cp'^ sin 29%) , (32) 

where 9% is the mixing angle in matter at the surface of the Earth, 9% = 
9"^{xf), Xf is the surface coordinate, and 



dx 



d9'^2ix 



dx 



sind'^'^ix). 



(33) 



Here x is the central point of the trajectory, and the expressions for probabili- 
ties are valid for a symmetric (with respect to x) density profile. The adiabatic 
half-phase is given by 



dr^ 



cos 2^12 — 



2EV{r)c 



-13 



+ sin^ 2^12 (34) 



21 



and in eq.(31, 32) 0"'^ = (f)'"^{xf). 

Using the amplitudes (31, 32) we obtain the following expressions for P2, -R2, -^2: 



P2 = cos^ I0 sin^ 29% sin^ 0"'^ + - sin 27^ sin 49% sin (p"^ + sin^ Ig cos^ 29%, 



Ri = -- 



cos^ 1q sin 49% sin=^ 0'*'^ + sin 2/^ cos 49% sin 0"'^ - sin^ 1q sin 40 



cos^ le sin 29% sin 20"'^ + sin 2Iq cos (pf^ cos 29% 



(35) 
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These expressions are reduced to the expressions in (30) with substitutions 
9^ and ^"'^ — * 0™ if = which would correspond to a perfect adia- 

baticity. From (35) we find formulas with the first non-adiabatic corrections 



P2 ~ sin^ 2^°2 sin^ <p'"^ + h sin 4^°2 sin 

i?2 ~ - - sin 4^°2 sin^ 0"'^ - h cos 4^°2 sin 0''^ 

72 f« ^ sin 2^°2 sin 20"'' + 7^ cos 20°2 cos 4)"^. (36) 

For completeness in the Appendix A we present also results in the second 
order of the Magnus expansion which provides better than 1% accuracy for 
all neutrino energies. 



The results of this section can be used immediately for high energy neutrinos 
E > 150 MeV and for trajectories far from horizon, when oscillations in the 
atmosphere can be neglected. 



4 Oscillation effects in the atmosphere and inside the Earth 



Above horizon the oscillations occur in vacuum (we neglect density of the 
atmosphere). The evolution of neutrinos can be considered again in the z/' 
basis, since the connecting matrix U' depends on the vacuum mixing angles 
only. 

Again, in the u' basis the third mass state decouples and for two other states 
the 2x2 block of the S-matrix is given by 



(2) I c,^ + « cos 2^125,^ -i sin 2^125^ ^ 
—i sin 2^125^ — i cos 2^125,^^ 



(37) 



Here = sin0, = cos0 and the half-phase is defined in (7). The oscil- 
lation probabilities in the flavor basis are given by the expressions (27 - 29) 
with P2 — sin^ 2^12 sin^ (p, R2 — sin 4^i2 sin^ 0, and 72 = | sin 2^12 sin 20. 

For trajectories below the horizon one needs to take into account oscillations 
in two layers: the atmosphere and the Earth. For these trajectories the total 
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S— matrix in the i/'— basis equals 



■>tot 























^33 



(38) 



where S' is given in (14) and Sa describes evolution in atmosphere. 

Using the expressions for the ampUtudes (85) from the Appendix B and prop- 
erties of the amplitudes for a single symmetric layer (17) we obtain the rela- 
tions 

(39) 



A —A* 



-A 



The flavor oscihation probabihties can be computed according to P(t'a 
Up) = \ {U'StotU'^)i3a\'^- As a result, we obtain formulas similar to those in (27 
- 29) for a single layer with, essentially, substitution A'^^ — ^ 



P{^e ^e) — Ci3|^ee|^ + *13) 

P{j^H T^e) = -Si^S2ze~'^Aee + C23A 



-13 



-Sl3S23e* Aee + C2^A 



fie 



I 2 2 2 
+ *13'^13*23' 

I 2 2 2 
+ '^13'^13'^23' 



(40) 
(41) 

(42) 



cla^nn - ^^13 sin 26'23lm(e*'^^e/,) + sl^sl^A^e + ctasta-i^^) 
In the last equation we used the properties (39). 

Since the two-layer density profile is not symmetric, we have inequality A^^ ^ 
Ae^, and furthermore, these amplitudes are not pure imaginary. Consequently, 
the probabilities (27 -43) are expressed now in terms of 5 different functions 

and not three, D, as in the case of symmetric profile. We introduce them as 



P — \^tie\ — l^e^til ) Ren — ^^{.-^eix^ee) i -^e/i — ■'■ni(v4.g^Aee) , 

i?^e = Re(i;,iee) , he = Im(i;,iee) ■ (45) 

From unitarity of the S— matrix we have = l^eeP = 1—P and A^eA*^^ = 

— Ag^A*^. In turn, P, Re^, /g^, R^e and J^e are the functions of P2, R2, h and 
vacuum the oscillation parameters (p and 6*12. We present these expressions in 
the Appendix B. 
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In terms of D— functions (45) the expressions for probabilities arc similar to 
those for 1 layer (27 - 29) with, however, certain differences. Using (40 - 43) 
and definitions (45) we obtain 

P(i/e-^e)=4(l-P)+4, (46) 
Pil'u, l^e) = ci^cl^P - S13C13 sin 26*23 (cos 5i?eM + single/,) + Si3Ci3S23(2 - P). 

(47) 

Then 

P{^e ^n) = C13C23-P - Si3Ci3Sin26l23(cOs5i?/,e " Sin 5/^e) + Sl3C?3«23(2 " ^)- 

(48) 

The most significant change is in the — > probability: 



P{y^. -^iy^) = l- 0.5 sin^ 2^23 - + 

+ Si3sin26'23(c23 - sl^sl^)[cos6{Ref, + R^,e) + sm6{Ie^ - I^,e)] + 

p { 1 

+ S?3 Sin^ 2^23^ |1 + p(^-Y_p^ [C°^ 2(5(i?e^i?^e + + 
+ Sin25(i?^e^^ — Retinue)] + (-^e//-^/ie " -Re/i-R/ie)| — 

- 244 + 44(2 -P). (49) 
In the limit (no oscillations in atmosphere) we have 

Pen — P^e — P2, = Ifxe = P = P2, (50) 

and the formulas (46 - 49) are reduced to the one-layer formulas of the pre- 
vious section. Notice that functional dependence of the probabilities on the 
parameters S13, 5 and ^23 is practically the same as in the one layer case. 

Let us consider expressions for i?e^, R^ie, hfi, I^e and P in the constant density 
approximation (see Appendix C.) They can be presented in the following form 

D,^ = J^UcT + 0) + sin Xi{0T2, O12), (51) 

where = O^ — Ou- So, D? can be written as the corresponding functions for 

one layer of matter with total phase (f)'^+(f) plus corrections related to difference 
of mixing angles in matter and vacuum. Since at low energies A^ <^ 612, the 
second term in (51) can be considered as small correction. In the first order 
in A^ the expressions (95) become 
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P" ^ sin^ 2^5^ sin\(f)"' + 0) - 2Ae sin 4^i2 sm(0"* + 0) sin cos 0"^, 
^e/. ~ - ^ sin 4^™ sin2(0'" + (j)) + AO [- sin^ 2^^^ sin 20 sin 20"^+ 



+ 2 sin^ 0(cos4^i2 cos^ 0™ + sin^ 0"^)] , 

~ - ^ sin 4^^ sin2(0"^ + 0) + [cos^ 2^^^ sin 20 sin 20"*+ 



+ 2 sin^ 0(cos4ei2 cos^ 0"* - sin^ 0™)J , 
11^ ^ ^ sin 2^™ sin 2(0™ + 0) - cos 29'^^ sin 20, 
^ sin 2^™ sin 2(0™ + 0) + 



+ A^ cos 2^^ - sin 20 cos 20™ + 2 sin^ sin 20^ 



(52) 



At energies £■ < 60 MeV the magnitude of second term does not exceed 0.1. 
For high energies the difference A^ becomes large. However in this case the 
oscillation length increases and becomes smaller, so that the oscillations in 
atmosphere can be neglected. 

Summarizing, final expressions for the flavor probabilities are given in eqs. 

(46), (47), (48) (49) with P, Re^, Je^, i?^e, /^^e as functions of P2, R2 and h, 
and the vacuum oscillation parameters presented in the Appendix B. In turn, 
the precise semi- analytical expressions for P2, R2 and I2 are given in (35) for 
the first order Magnus expansion and in the Appendix A - for the second 
order. 



5 Averaging of oscillations and integration over the zenith angle 

Results for the oscillation probabilities obtained in the previous section are 
simplified substantially after integration over the zenith angle. In what follows 
we will consider detection of neutrinos by the charged current interactions with 
nucleons, where the information about neutrino direction is essentially lost. 
Consequently, observables (at least in the first approximation) are given by 
integration over the zenith angle. To a good approximation the neutrino flux at 
low energies does not depend on the zenith angle Q^, and therefore integration 
over the angle is reduced to averaging of probabilities over 0,^. 

In the constant density approximation the probabilities averaged over the 
oscillation phase equal 



where subscript means averaging over the phase. For antineutrinos one should 




(53) 
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substitute — >■ As a consequence of /| = the CP-odd terms of prob- 
abilities (which lead to the CP-asynimetries) become zero, as is expected, 
since this would correspond to averaging over all oscillation phases (we have 
already averaged over the 1-3 phase). The real part i?2 determines the first 
order correction due to the 1-3 mixing. 

For very low energies we have 6^ ^ 9i2, and consequently, (-R2)</> = {^2)1^ = 
— I sin 4^12 ~ —0.19. For antineutrinos with increase of energy the angle 
9'^2 decreases. As a result, (i?2)(/. < in whole energy range. The absolute 
value |(-R2)</)| first increases, reaches maximal value 0.25 when — 22.5° (at 
E ~ 150 MeV) and then decreases. In contrast, the angle increases, and 
therefore (-R2)</> first decreases, becomes zero in the resonance {E ~ 100 MeV) 
and then changes the sign: (_R2) > 0. It increases with energy till E ~ 300 
MeV and then decreases again. 

For the two layer case we obtain after averaging over the phases 



{P")^ = \ sin^ 2^™ - 2^™ - sin^ 2^12) + +\ sin^ 2(^12 - ^^2), 

(^^0 —\ 4^2 + \ [sm4er2 - sin 4^12 + sin4(^- - ^12)] , 

{RU)4> = -\ sin4^^r2 + \ [sin 4^™ - sin4^i2 - sin4(^™ - ^12)] , 
(4y^=U^e)0 = O. (54) 



At very low energies, when the difference of mixing angles in vacuum and in 
matter is very small, we find from (54): 



(P)^ R:ilsin22^r2-^A^sin4^i2, 
(i?e/.)^ = -\ sin4^^"2 + A^cos' 2^12, 

(i?^e)<A = -\ sin4^r2 - A^sin^ 2^12. (55) 



These averaged functions (A)</. determine the number of low energy events. 
Indeed, to find the number of events we need to integrate the probabilities 
folded with the neutrino fiuxes and cross-sections over the zenith angle 
(the angular variables, in general). In turn, the oscillation probabilities are 
hnear functions of Dj. Therefore we deal here with integrals 

J, = /(i(cose,)A(cose,)F(e,). (56) 
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At low energies the fluxes do not depend on Q^, and consequently, 

Ji = 2F (A(cose,)). (57) 

Here (...) without subscript denotes averaging over the zenith angle Oj^. The 
functions Di depend on cos Gj, via the oscillation phases 4> and (pm- Apparently 
(f) oc L^(cos6,^), and in the constant density approximation, 0^ oc oc 
cos©,^. Therefore averaging over cosOj/ is equivalent to averaging over the 
phase (f)jn as we did above. The situation is different for the terms which depend 
on the vacuum phase. Indeed, according to (4) the dependence of on cos 9;^ 
is non-linear. In particular, for trajectories not very close to the horizon oc 
1/cosO,^ and therefore integration over cos 9,^ is not reduced to averaging 
over the phase. For estimations, the effect of oscillations in atmosphere can be 
taken into account by the additional factor k as 

(D,(cos9,)) « ^ 1 (i(cos9,) D, = |(D,)<^, (58) 

where (Di)^ refers to averaging over the phase. The factor k equals 1, if os- 
cillations in atmosphere can be neglected. For low energy bins (£■ < 15 MeV) 
K, can be as big as 1.2; for E — 20 MeV we have k ~ 1.1 and k, approaches 1 
with increase of energy. 



6 Neutrino fluxes at the detector 



The z/g— flux at the detector with oscillations taken into account, F^, can be 

written as 

Fe = F^P{u, ^ v,)+FlP{v^ ^ V,) = F° [P{v, ^ u,) + rP{u^ ^ v,)] , (59) 

where r is deflned in (1). For antineutrinos we have similar expression with 
substitutions P^P,F^F,r^f. 

Inserting the probabilities Pee and P^e from (46) and (47) into (59) we obtain 
the expression for the relative change of the i/g— flux due to oscillations: 



F 

-^ = 1 + (rc23 - 1)P - rsi3Ci3 sin 26'23(cos 5 Reu, + sin 5 hfj) 



-2s 
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1 - r4) + P(r - 2] 



+ 4(l-r4)(2-P). 



(60) 



It coincides (up to the sign of 5 and the corresponding change of D ^ D) with 
the expression in our previous paper [15]. Notice that the corrections of the 
order are suppressed by the "screening" factors which are zero for maximal 



2-3 mixing and r — 2. There is no corrections of the order s 



13- 
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The fy^— flux at the detector with osciUations taken into account, F^, can be 
written as 



Then the ratio of fluxes with and without oscillations equals 



(61) 



1 - i sin^ 2^23 - C23-P \cl^ - 



-13 



--23 



Si3 sin 2^^23 { cos 5 



R 



'/le 



— sin 5 



^ + 4(^e. + m+0(«?3). 



(62) 



Averaging over the zenith angle vanishes the imaginary parts. Notice that the 
linear in S13 term is proportional to i?e/i + Rue ~ 2i?2- 



7 Electron neutrino (antineutrino) events 



We will discuss here (mostly for illustration) a detection of the low energy 
neutrinos in the water Cherenkov detectors. In these detectors the electron 
neutrinos are detected via the quasi-elastic scattering 



i>e + p ^ e+ + n, Ve O ^ e+ + AT, V° O ^ e' + F. 
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(63) 



The energy of positron (electron) is practically uniquely related to the neutrino 
energy: Ee = - 1.293 MeV, (for ^H), E^ = E^ - 15 MeV (for ^^O). The 
number of e"*" events (63) has been computed as 

A^eoc ^ I dE,dE,d{cose,) F,{E,,e,)^-^{ee,e,,E,)e{E,), (64) 

i=p,0 6 

where Fg is the atmospheric Pg-flux at the detector given in (60) for 9^; the 
fluxes Fg and without oscillations are taken from Rcfs. [3,4,6]; dai/dEe 
are the differential cross-sections [17], e{Ee) is the detection efficiency. \E' is 
the "dispersion" function which describes deviation of the lepton zenith angle 
from the neutrino zenith angle (for details see Ref. [18]). Summation proceeds 
over scattering on Hydrogen and Oxygen. For we use similar expression for 
number of events but consider the reaction on Oxygen only. 

Results of computations of the energy spectra of u^— and i/g— events for 
different values of 2-3 mixing are shown in fig. 4. Qualitative and to a large 
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extend quantitative understanding of these results can be obtained from the 
following semi-analytical consideration. The number of events can be presented 
as 

^e=(^)iVr, (65) 

where is the number of events without oscillations; the ratio of fluxes Fg/ 
is given in (60) and averaging proceeds over the neutrino energy bin and the 
zenith angle. 

The relative change of the spectrum can be written in the following way: 

§ = l + e(°)+eW+e(^) + ..., (66) 

e 

where 

4°^ = (PiE)) [(r(E))4 - l] (67) 
is the correction due to the oscillations driven by the 1-2 mixing and split for 
Si3 = 0; 

ei'^ = -sn cos 5 sin 2^23 (r) {Re^) (68) 

is the first hnear correction due to 1-3 mixing, etc.. Here {P{E)), {Re^{E)) 
and {r{E)) are the probabilities and the ratio of fluxes averaged over the 
zenith angle and energy bin correspondingly. We take (ri?2) ~ (r)(i?2)- For 
antineutrinos P2 — > ^(-E) and r ^ f. 

In fig. 4 (left panel), we show distortion of the energy spectrum of positrons 
produced by the atmospheric for S13 = 0. All the features of the curves 
can be traced from eqs. (66, 67). The averaged probability {P2{E)) slightly 
decreases with energy. For estimations one can use the oscillation probability 
for a layer of constant density, so that {P2{E)) ~ 0.25k sin^ 2^^, where k takes 
into account oscillations in atmosphere (58). 

For the antineutrino channel sin^ 29^ decreases with increase of energy. Esti- 
mations give {P2{E)) = 0.089, 0.038, 0.013 for E = 20, 40, 80 MeV correspond- 
ingly. For a given cig the change of eg°' reflects this decrease of {P2{E)) and 
change of {f(E)). In general, ef^ decreases with energy. Distortion of spec- 
trum can be characterized by Aeg*^^ = ee(100 MeV) — ee(10 MeV). According 
to the figure the distortion Aeg°^ is stronger for small C23, it decreases from 
~ (8 - 9)% for c^3 = 0.35, down to ~ (2 - 3)% for 4^ = 0.65. Thus, if the 
shape of spectrum can be predicted with accuracy ~ 1%, one can measure the 
deviation of 2-3 mixing from maximal by studying distortion of the spectrum 
in a way which does not depend on uncertainties of the absolute value of neu- 
trino flux. One possibility is to compare the integrated signal below and above 
53 MeV. 

For cig > 0.6 the dependence is not monotonous: with decrease of eg first 
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Fig. 4. Left panel: Distortion of the positron energy spectrum due to oscillations of 
the atmospheric Pg- Shown is the ratio N^/N^ as function of the positron energy for 
different values of cos^ 023- The Bartol fluxes and values of oscillation parameters 
sin^ 26*12 = 0.82, Am|i = 7.3 x 10"^ eV^ and S13 = used for calculations of the 
histograms. Both interactions with and ^^O in a water Cherenkov detector are 
taken into account. Right panel: The same as in the left panel for the atmospheric i/g 
and the electron spectrum. The interactions with ^^O in water Cherenkov detectors 
are taken into account. 

increases and then decreases. For C23 > 0.6 the oscillation effect is positive 
eg > in whole energy range, whereas for < 0.42, ef^ < everywhere, and 
for C23 ~ 0.5 the correction changes the sign. 

With change of C23, Aef^ = ee(0.65) — eg(0.35) = (11 — 12)% for large energies 
{E ~ 90 MeV), and Aeg = 16% in the low energy bin. So, if the absolute 
flux is known with accuracy few %, the 2-3 mixing can be determined by 
just measuring the total number of events. Variation of the signal with C23 
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Fig. 5. Distortions of the energy spectra of events for the neutrino fluxes from 
different computations. The full and dashed lines correspond to two different values 
of cos'^ 023- Left panel: Ue, Right panel: Ve- 

is stronger here than in the sub-GeV sample {E > 100 MeV), although the 
number of events is smaller. 

In fig. 4 (right panel), we show the distortion of spectrum of the e— like events 
induced by scattering of z/g on ^^O. For neutrinos the flavor ratio r is sys- 
tematically smaller than for antineutrinos. Therefore the histograms shift to 
smaller values of ratio Ng/N^. The suppression of signal due to oscillations 
is stronger here. A character of the distortion is rather similar to that in the 
z/— case, and again, it can be traced from eqs. (66, 67). Now {P2{E)) changes 
with energy weaker. Indeed, k decreases, whereas sin^ 29m increases with the 
energy increase and the two changes partly compensate each other. For a given 
we obtain Ae(°) = e(o)(100 MeV) - e(°)(10 MeV) ~ 6% for all values of 
in the interval 0.35 — 0.65. For fixed energy, with change of C23 one has 
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Fig. 6. Effect of the 1-3 mixing on the spectra of events induced by the electron neu- 
trinos and antineutrinos. Shown are relative spectra of the e— like events for different 
values of 1-3 mixing. We take cj^ = 0.5, sin^ 2(9i2 = 0.82, and Am|^ = 7.3 x 10"^ 
eV^. Left panel: Ug, right panel: v^.. 

e(°)(0.35) - eW(0.65) ~ 15% for all energies. 

In fig. 5 we compare the distortions of spectra of events computed with the 
neutrino fluxes published by different authors. This quantifies the present 
theoretical uncertainties. Fluxes from [4] (Honda) and [3] (Bartol) lead to 
rather similar distortion although according to [4] riE) is flatter than in [3] 
in the range > 60 MeV, but below that energy r{E) decreases sharper (see 
fig. 2). For large C23 this leads to non-monotonous change of eg with maximum 
at = (60 — 70) MeV, as we have marked before. 

Corrections due to 1-3 mixing described by e^^-* + e^^-* + ... are shown in fig. 6 
for antineutrinos (left panel) and for neutrinos (right panel). All the features 
of the figures can be immediately understood using the analytic expressions 
derived above. Recall that due to integration over the zenith angle J2 becomes 
negligible (see eq. (53)). Moreover, as we marked before, the corrections of 
the order are additionally suppressed. Therefore in the first approximation 
the corrections are given by e^}^ in eq. (68). This term weakly depends on ^23- 
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{Reij) can be estimated taking the constant density approximation: 

(i^e^)- W~-|sin4^-- (69) 

For directions below the horizon the corrections due to oscillations in atmo- 
sphere are additionally suppressed by cos^2^i2 ~ 0.18 (55). 

Consider first the effect for antincutrinos. In the lowest energy bin we have 
(f) = 1.7, n ~ 1.2, and {Ren) = —0.085^. Therefore the linear in S13 correction 
equals 

ey'(lOMeV)^ 0.021 (70) 

With increase of energy both (^2) and (f) increase and at £■ = 100 MeV 
(0.5 sin 4^^ ^ 0.48): 

4^)(100 MeV) 0.033(^1^) (71) 

in agreement with the results in fig. 6. The corrections are CP-even being of 
the same sign for neutrinos and antincutrinos. 

For neutrinos we have (r) = 1.66 in the lowest energy bin, and consequently, 

e«(10 MeV) 0.025(^1^). (72) 

The correction decreases with energy due to decrease of (Ren)- For £■ = 60 
MeV we have (r) = 2.0, (Re/,) = -0.035, and therefore 

eW(60 MeV) 0.009(^1^). (73) 

In the resonance i?e^ ^ 0. Therefore, in the region around ~ 90 MeV the 
corrections due to 1-3 mixing are suppressed additionally, as can be seen from 
fig- 6. 

Comparing results from different energy regions we can identify the 1-3 mixing 
effect (induced interference). Notice that corrections due to 1-3 mixing weakly 
depend on the 2-3 mixing which allows us to disentangle the effect of 1-3 
mixing, and the effect due to the 1-2 mixing which is proportional to the 
deviation of 2-3 mixing from maximal. 

Corrections proportional to sl^, 

6(^)~24[((r)4-l) + (P)(2-(r))], (74) 
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are small: < 0.005 for 5^3 = 0.016. Thus, the degeneracy of sf^ and C23 param- 
eters can in principle be resolved using the energy dependence of the effect for 
neutrinos, however for antineutrinos the energy dependence is weak. 



8 Muon neutrino detection. Invisible muon decays 



Muon neutrinos can be observed via the decay of invisible muons. Muons with 
energies below 160 MeV do not produce signal in water Cherenkov detectors. 
They quickly lose energy, stop and decay at rest emitting electron (positron). 
In turn, most of these muons are produced in the detector in the quasi-elastic 
interactions + N ^ N' + fi. The effective (reconstructed) energy spectrum 
of these neutrinos has maximum at 0.16 GeV and extends to 0.25 GeV [8]. 
In this energy range ki Fp, Fg > F^ and the flavor ratios equal r = 2.03, 
and f = 2.23 [6]. Typical energies of neutrinos are high enough so that for 
estimations we can neglect the oscillations in atmosphere and use functions 
D. 

Consider effects of oscillations on the number of electrons/positrons from in- 
visible muon decay. Apparently, the shape of energy spectrum of events is not 
influenced by oscillations and it has the standard form with maximum at 45 
MeV. In flg. 7 we show the ratio of number of the invisible muon decays with 
and without oscillations, N^/N^, as function of C23 for different values of S13. 
All the features of this flgure can be traced from the following semi-analytic 
consideration. According to (62) we have approximately 

^^l-lsin^2^23 + i°)+6«+e(') + ... (75) 
fj. 

The zero order (in S13) correction equals 

- -(p^yis (4 - ^) • (76) 

The linear in S13 corrections is given by 

e'j}^ « sia sin 2^23 cos S {R^) (^2^ - . (77) 
The quadratic in S13 corrections can be written as 



e?^ ^ s^3 <! sin^ 2^. 
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1 //l 



-24 + iy(24-(P2))} 



(78) 
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Fig. 7. Oscillation effects on the total number of invisible ;U— decay events induced 
by antineutrinos {left panel) and neutrinos (right panel). Shown is the dependence 
of N^/N^ on C23 for different values of 1-3 mixing. Dashed line corresponds to the 
vacuum i/^ — Vr oscillations (zero 1-2 and 1-3 mixings). 

In the constant density approximation we have 

11^3 = 1 - (1 + cos' 2^™) sin' 0"^, (79) 
and averaging over the zenith angle gives 

(^)4(' + ^-^^^s). (80) 



The first term in (75) is just the averaged vacuum z/^ — Ut oscillation prob- 
ability. This probability (black curve) is symmetric with respect to maximal 
mixing, C23 = 0.5. As expected, at this value of C23 one has maximal suppres- 
sion, 0.5. 

Corrections due to the oscillations driven by the 1-2 mixing, e^^\ are asymmet- 
ric with respect to 0^3 = 0.5. The corrections are zero at C23 = 1/ (r). They are 
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positive at smaller 0^3 and negative at C23 > 1/ (r). Maximal relative effect is 
~ 5% suppression at C23 = 0.65. The effect is weaker (< 3%) for antineutrinos 
which is related to smaller probability P2 and larger value of f only partly 
compensate the diflference. 

The first order correction due to the 1-3 mixing, e^^^^ is given in (77). For the 
average energy of neutrinos producing invisible muons, E = 160 MeV, we have 
{R2) = 0.04 and (^2) = —0.125. Notice that the energies of neutrinos which 
generate the invisible muons are above the resonance, and therefore R2 has 
opposite to R2 sign. Then for c^^ — 0.65 we obtain 

For antineutrinos the corrections are much larger and of opposite sign: 

^ -0.0127 ( '-^] . (82) 
" \ 0.126 J ^ ' 

This is in agreement with the results of fig. 7. 

In contrast to the direct - channel, here 5^3 corrections are not suppressed 
by the screening factors and turn out to be of the same order as the linear 
corrections for not too small 5^3. Using (78) we get for the quadratic term and 
cos (5 = 1 

e(f) ^ 0.0085 ( , (83) 

" Vo.oi6y' ^ ^ 

and for antineutrinos the corresponding numerical coefficient equals 0.0054. 
The corrections are positive. 

Notice that the high energy neutrino fluxes responsible for the invisible muon 
production are affected by the solar activity weaker than fluxes at low energies. 
At £^ > 150 MeV the difference of fluxes during minimum and maximum of 
solar activity is (6 — 10)%, whereas at 30 MeV the difference is about 24 %. 
Variations at high energies are weaker. This can be used to disentangle the 
direct production and signal from the invisible muon decays. 



9 Future measurements 

High statistics of the sub-subGeV events can be obtained in future megaton- 
scale water Cherenkov detectors [19,20,21,22]. In these detectors the signals 
of four different neutrino fluxes considered in the previous sections can be 
identified in the following way: 
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1) . The Ve quasi-elastic scattering on protons can be detected performing tag- 
ging of neutrons in correlation to the positron detection. In SuperKamiokande 
this will be possible with Gadolinium in the GADZOOK's version [23]. 

2) . The T>n quasi-elastic scattering on protons has similar signatures: detection 
of neutron and positron from the muon decay. The difference from the previous 

case is that positrons appear in certain energy range and have known energy 
spectrum. So, the number of invisible decays can be extracted by fitting the 
energy spectrum. (Also neutrons will have different energy characteristics.) 

3) . The electron neutrinos v^s are detected by their scattering on Oxygen. 
Electron appearing in this process should not be accompanied by neutron. 
(In principle, one can perform tagging of nuclear transitions, e.g., detecting 
transition to exited state of nuclei.) 

4) . The u^— flux is detected via the invisible muon production on Oxygen. 
The standard decay spectrum of electrons and absence on neutron are the 
signatures of this interaction. 

Thus, the antineutrino-induced events can be disentangled from neutrino in- 
teractions by the neutron tagging. 

Certain conclusions can be drawn from comparison of different signals. The 
oscillation corrections to the Vf.— signal are substantially enhanced in compar- 
ison to the Pe corrections if < 0.5. Comparison of the a-nd Pg signals at 
> 60 MeV can reveal the effect of 1-3 mixing. 

Fig. 8 illustrates dependence of the spectra of the e— like events on the 2-3 
mixing for S13 = 0. The histograms show spectra of the e— like events produced 
directly by the electron neutrinos and antineutrinos (red and blue) , and from 
the invisible muon decay (black) for different values of cos ^23- We show the 
number of e-like events at HyperKamiokande (540 kton collected during 4 
years). The difference of numbers of events for C23 = 0.35 and = 0.65 
in each bin is about 2a (statistical error). Therefore using 9 bins one can 
distinguish two values of C23 at IScr level. 

The effect of 1-3 mixing is much stronger in the antineutrino signal than in 
the neutrino signal both for the electron and muon neutrinos (i.e., for both 
components shown in the fig. 9). The difference of histograms is essentially 
due to the linear sia-corrections. Qualitatively the size of the effect can be 
immediately inferred from eq. (68) for direct component and from eqs. (77, 
81, 82) for the invisible muon decays. 

For maximal allowed value 5^3 = 0.03 one would get 1.4 times larger difference 
of histograms than the one shown in fig. 9. 
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Fig. 8. The spectrum of neutrino and antineutrino events expected for 4 years of 
the HyperKamiokande (0.5 Mton) data-taking. The red and blue histograms are 
for the e-hke events directly produced by and i>e for two different values of 2-3 
mixing. Black histograms correspond to the e— like events from the invisible muon 
decays computed respectively for cos^ ^23 = 0.35, and cos^ 6*23 = 0.65 (dashed lines) 
and cos^ ^23 = 0.50 (solid line). 

Notice that with increase of 0^3 the number of invisible /i— decays increases 
for C23 > 0.5 and decreases for C23 < 0.5. The number of e-hke events always 
increases with C23. This can be used for measurements of C23. 

Tagging of neutrons allows one to detect total signal produced by antineutri- 
nos: Ve and v^. One can further disentangle the contributions from electron 
and muon antineutrinos by making the fit of spectral shape and using known 
spectrum of positrons from the muon decay at rest. This contribution is char- 
acterized by a single normalization parameter which can be then extracted 
from the data fit. 
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Fig. 9. Effect of the 1-3 mixing on the spectra of e— hke events. The total number of 
e— hke events induced by v^., Ve and invisible muons decays for the HyperKamiokande 
4 years of data taking is given as a function of the charged lepton energy. Shown 
are the histograms for 0^3 = 0.5, sfg = 0.016 and cos 5 = +1 and cos (5 = —1. 

If no neutron tagging occurs, one measures the total energy spectrum of elec- 
trons and positrons (see figs. 8 and 9). As follows from the figures, even in this 
case one can extract certain information on deviation of the 2-3 mixing from 
maximal and on the 1-3 mixing. 

As we have mentioned before, different components of spectrum have different 
dependence on the solar activity. This allows one, in principle, to disentangle 
the direct z/g— , i/g— signal, and the signal from the invisible muon decays (muon 
neutrinos) . 

High statistics of future experiments will open a possibility to play with other 
characteristics such as angular distributions, directionality, time dependence. 
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Also detection of the neutral current events can provide an important infor- 
mation. 

A very detailed study of the sub-subGeV events including determination of 
flavor and charge of the produced lepton will be possible with large volume 
hquid Ar detectors [24,25], [26]. 



10 Conclusion 

1. There are several new features which appear at low energies {E < 100 
MeV) in production, oscillations and detection of the atmospheric neutrinos. 
One of these features is substantial decrease of the flavor ratios with decrease 
of energy and their substantial deviation from 2 in very low energy bins. 
As far as oscillations are concerned, there is strong averaging of oscillations 
due to integration over the angular variables. The oscillations in atmosphere 
become important. For events induced by the electron neutrinos the energy of 
electron (positron) is practically related to the energy of neutrino. In contrast, 
for muon neutrinos detected via the decays of invisible muons, information on 
the neutrino energy is practically lost. These events are induced, mainly, by 
neutrinos in the energy interval (0.1 - 0.3) GeV. 

2. We performed both numerical and analytical study of the oscillations of 
sub-subGeV neutrinos. Dependence of the oscillation effects on sin 6*13, 5 and 
deviation of the 2-3 mixing from maximal are given explicitly. For the rest one 
can use either numerical results or semi-analytical formulas. 

3. For neutrino trajectories below the horizon one should take into account the 
oscillations in two layer medium: the atmosphere and the Earth. We have pre- 
sented the relevant oscillation probabilities and numbers of events as functions 
of sin ^13, S, 023 and in terms of P2, R2, h- For the latter we give precise semi- 
analytical formulas using the Magnus expansion as well as explicit expressions 
for the constant density. 

4. For the e— hke events at S13 = the corrections due to the oscillations 
driven by the 1-2 mixing can be as large as 10 — 15% depending on possible 
deviations of the 2-3 mixing from maximal. The oscillations lead to change of 
the total number of events and also to distortion of the energy spectrum. The 
distortion is stronger in the low energy bins. The effect (change of the slope 
of the spectrum in the interval 10 - 100 MeV) can be as large as (7 — 8)%. 

5. The 1-3 mixing leads to an additional, somehow smaller, effect. For maximal 
allowed values 5^3 ~ 0.03, it can reach ±4% at low energies, and ±6% at 
high energies. The energy dependence of the effect is stronger for Ug. The 
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linear corrections which dominate at low energies vanish in the resonance 
~ 100 MeV. For antincutrinos corrections shghtly decrease with energy. The 
corrections for antincutrinos are much larger than for neutrinos at £^ ~ 100 
MeV. 

6. In the energy interval (0.1 - 0.3) GeV the i/^-flux is detected via the decay of 
invisible muons. The total number of e— like events from these muons depends 
on the 2-3 mixing in the lowest order. For allowed values of the deviation of 
the 2-3 mixing from maximal the effect can be as large as 5%. Inclusion of the 

1- 2 mixing gives (for S13 = 0) maximum 5% effect for neutrinos and 3% effect 
for antincutrinos. This maximal value is realized at c^g = 0.65. For C23 < 0.5 
the effect is smaller. 

The 1-3 mixing produces comparable to 1-2 mixing corrections. The largest 
effect is for antineutrinos and c^g — 0.65: e^''^^ ~ 0.02. The linear in S13 cor- 
rection is strongly suppressed for neutrinos. The quadratic corrections are not 
neghgible and can be of the order of 0.01. 

7. There is strong degeneracy of parameters S13 cos 5 and c^g, especially for the 
invisible decays. Inclusion of the 1-3 mixing effects only slightly enlarge the 
region of possible values of ratio Nn/N^. 

8. The signals from interactions of 4 different fluxes u^^ ^'^ '^ei ^'^ can be disen- 
tangled by tagging the accompanying neutrons, studying the shape of energy 
spectrum as well as the time dependence of signals. Confronting events of dif- 
ferent types as well as events at low and high energies allows one to reduce 
the degeneracy of oscillation parameters. 

9. The number of presently detected events (at SuperKamiokande I) due to 
the interactions of very low energy neutrinos is about few hundreds which 
provides an accuracy of determination of the oscillation parameters not better 
than 10%, and an additional problem is uncertainties of neutrino fluxes. So, to 
study oscillation effects discussed in this paper one needs much larger statistics 
which can be achieved with the Megaton-scale detector. 

10. The low energy atmospheric neutrinos can be used to measure deviation of 

2- 3 mixing from maximal, the 1-3 mixing and the phase 5. They can be used 
to search for new physics. Apparently for the latter knowledge of the standard 
oscillation effects computed in this paper is necessary. Understanding of fluxes 
of these neutrinos is also important for future studies of the relic supernova 
neutrinos. 



31 



Acknowledgments 



The authors are grateful to G. Battistoni and T. Stanev for discussions of 
the atmospheric neutrino fluxes at low energies. O. L.G. P. thanks FAPESP, 
CNPQ and ICTP for support. 



Appendix A: Probabilities P2, R2 and I2 in the second order of the 
Magnus expansion. 

The probability functions D in the second order Magnus expansion can be 
found using definitions (24) and results of [16]: 



cos Itot sin (f)""^ - ^ sin hot cos 

^tot 

cos hot sin 0" - — sin hot cos 

^tot 



P2^ 



i?2 = -^sin4C 



."'^'j sin2^^ + ^sin/i„tCos2^^ 

/ ^tot 



\2 



— sm hot 

^tot^ 



- COS 46*^ cos hot sin 0"" - — sin hot cos 0" j — sin hot, 

V ^tnt ' hot 



= \ cos^ hot sin 2^^ sin 20"'* + sin Ihot cos 2^^ cos (f^ + 

2 21 tot 

sin 2^ cos - ^ sm' hot sin 2^ sin 20"'*+ 



^^cos2Csin0'^'* 



-'tot 



Here 



and 



hot = \/-^| + lee^ 



he = 4 dx /_ dy 

Jx Jx 



d9m{x) 
dx 



d9m{y) 
dy 



sin0'^'*(a;) cos 0'^'*(7/). 



(84) 



Appendix B: Functions D for two layers. 



According to eqs. (38, 14, 37) the elements of ^'-matrix in the two layer case 
equal 



32 



Aee = ^ee(c0 + 1 COS 261284,) - lA'^^j^ sin 26li2S0, 

A^^ = ^e^(c0 - i COS 2^125^) - iAgg sin 2^i2S<^, 

= ^e/.(c0 + i COS 2^i2S0) - lA'^^ sin 2^125,^, 

^^^ = A^^(c<^ - icos26'i2S0) - iAg^sin26'i2S0, (85) 

and A33 = ea;p[— i2(0?^ + 03)], where (j)^ is given in (15) and 03 is the phase 
acquired in atmosphere. We present here the probabihty functions D defined 
as 

(see also (45)). Using expressions for the amphtudes (85) we obtain 

P = P2 + R2 sin 46*12 sin^ + /2 sin 26*12 sin 20 + (1 - 2P2) sin^ 26*12 sin^ 0, 

Ren = -^2(1 - 2 cos^ 26*12 sin^ 0) - I2 cos 26*12 sin 20 - Q - P2^ sin46*i2 sin^ 0, 

lef, = 72(1-2 sin^ 0) + R2 cos 2^12 sin 20 + Q - P2) sin 2^12 sin 20, 

R I f 1 R^\ 

R^e = R2{l-2 sin^ 2^12 sin^ 0) H — |^ sin 2^12 sin 20 - - - ^ sin 4^12 sin^ 0, 



//.e = /2(1 - 2 sin^ 2^12 sin^ 0) + sin 4^12 sin^ + - - ^ sin 2^12 sin 20. 

P2 \^ P2 J 

(87) 

The hmit P2 = R2 = I2 = ^ corresponds to oscillations in the atmosphere 
only (trajectories above the horizon). In this case the probabilities are given 
by the last terms in each expression which agree with one layer results (30). 



Appendix C: Constant density case 



In the constant density approximation the evolution matrix in matter is given 
by 

c^ + icos2e^s^ -isin2e^s^ ^ 



S' 



-i sin 29VJ,sT 



- i cos 2e?ksT 



where = sin0"*, = cos0"*, and 0"* is the oscillation phase in matter. 
The total evolution matrix in two layers (vacuum and matter) is 5"*°* = S'Sa- 
Then in the basis v' the amphtudes for two layers, A^p — [S^°*]ap, equal 



iee = C^C^ - COS 2(^^2 " ^12)s^S^ + «(cOS 2^i2S^C^ + COS 2e^^S^C4,). 

A,^ = sin 2(^12 - ^r2)«^«<^ - ^(sin 2^125^0^ + sin 2^-^-0^), (89) 

A —A* A — _ 4* 
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The functions D (86) can be presented in the following form 

= sin^ 2^^ sin^ 0"* cos^ + sin^ 2^12 sin^ cos^ 0"^ + (90) 
+ ^ sin 2^^^ sin 2^12 sin 20"* sin 20 + sin^ 2(^™ - 612) sin^ 0"* sin^ 0, 



2Rl^ = - sin 4^™ sin^ 0™ cos^ - sin 4^12 sin^ cos^ 0™ - 

- sin 2^™ cos 2^12 sin 20"* sin 20 + sin 4(^^ - 612) sin^ (jT sin^ (91) 



27^^^ = sin 2^™ sin 20™ cos^ + sin 2^12 sin 20 cos^ 0™ - 

- sin(4e;^ - 2^12) sin 20 sin^ 0"* - sin 2^^ sin^ sin 20"*, (92) 



2R1^ = - sin 4^^ sin^ 0"* cos^ - sin 4^12 sin^ cos^ 0"* - 

- sin 2^12 cos 2^^^ sin 20"* sin 20 - sin 4(^^ - 612) sin^ 0"* sin^ 0, (93) 



27^^ = sin 2^^ sin 20"* cos^ + sin 2^12 sin 20 cos^ 0"* - 

- sin(4ei2 - 2^^^) sin^ sin 20"* - sin 2^12 sin 20 sin^ 0"*. (94) 

Notice that in the limit — > the first term in each expression corresponds to 
the probability of oscillations in matter, the second term gives the probability 
in vacuum (0^ 0). These two terms are related by the interchange of 
the vacuum and matter characteristics: 6m ^ 0, (j)m ^ 0- The rest is the 
interference of the oscillation effects in the two layers. The first two terms are 
the same in and 7?e^ as well as in 7^e ^-nd 7e^; in turn, these probabilities 
differ by the interference terms. The probability P is symmetric with respect 
to interchange of the matter and vacuum characteristics. 

For our discussion it is convenient to write the probabilities in the form given 
in eq. (51): 



34 



= sin^ 2^^ sin2(0™ + 0) - (sin^ 2^^ - sin^ 2^12) sin^ cos^ 0"" - 
- ^ sin 2^™ (sin 2^^ - sin 2^12) sin 20 sin 20"" + sin2(^™ - Ou) sin^ sin^ , 

R'e^ = -^ sin 4^^ sin2(0"^ + ~ ^ 2^™ (cos 2^12 - cos 2^^^) sin 20 sin 20"" + 
1 



+ -sin^ 



(sin 4^^2 - sin 4^12) cos^ 0™ + sin 4(^™ - ^12) sin^ , 



R^^^ = -- sin 4^™ sin^ (0"^ + 0) - 1 cos 2^^ (sin 2^12 - sin 2^^^ ) sin 20 sin 20™ + 

1 



+ -sin^ 



(sin4^;"2 -sin4^i 



,2 jm 



iu) cos 

= ^ sin 2^™ sin 2(0"" + 0) + ^ sin 20 
+ (sin 2^^2 - sin(4^™ - 2^12)) sin^ 



- sin 4 (^;^ -^12) sin 



2 



(sin2^i2 -sin2^^)cos2 0""+ 



m 



= ^ sin 2^™2 sin 2(0"" + 0) + ^(sin 2^12 - sin 2^^2) sin 20 cos 20"" + 
1 



+ ^ [sin 2^5^ - sin(4^i2 - 2^™)] sin2 0sin20"". 



(95) 
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